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We report on an experimental demonstration of laser wakefield electron acceleration using a 
sub-TW power laser by tightly focusing 30-fs laser pulses with only 8mJ pulse energy on a 100 
//m scale gas target. The experiments are carried out at an unprecedented 0.5 kHz repetition 
rate, allowing "real time" optimization of accelerator parameters. Well-coUimated and stable 
electron beams with a quasi-monoenergetic peak in excess of 100 keV are measured. Particle-in-cell 
simulations show excellent agreement with the experimental results and suggest an acceleration 
mechanism based on electron trapping on the density downramp, due to the time varying phase 
velocity of the plasma waves. 

PACS numbers: 41.75.Jv, 52.38.Kd, 52.65.Rr, 29.25.Bx 



Since the concept of laser driven plasma accelerators 
was first proposed by Tajima and Dawson [1], advances 
in high-power ultrafast laser technology have enabled 
successful production of energetic electron beams in nu- 
merous wakefield acceleration experiments[2 -4 . Plasma- 
based particle acceleration holds significant promise for 
future compact sources of relativistic electron beams be- 
cause of the large acceleration gradients plasma can sus- 
tain relative to conventional radio frequency cavities. A 
high intensity laser pulse propagating in an underdense 
plasma generates large amplitude plasma waves with 
phase velocities close to the speed of light. Under certain 
conditions, electrons can be trapped in the waves and 
accelerated to relativistic energies. Recent progress has 
demonstrated that using ultrashort laser pulses, relativis- 
tic electron beams up to GeV energies with low energy 
spread and small divergence can be produced in a highly 
nonlinear yet stable mechanism known as the "bubble" 
regime [5]. 

Accelerating electrons in the "bubble" regime requires 
a laser pulse that is both intense (with ao > 1, where 
ao = eA/rrieC is the normalized vector potential) and 
short (with pulse duration r < 27rc/ujp^ where ujp = 
Y^e^ne/meCo is the plasma frequency). In typical ex- 
periments, the intense laser pulse is focused onto the 
edge of a supersonic gas jet with a matched spot size. 
The densities of such gaseous targets indicate that ul- 
trafast laser systems with pulse energies on the order of 
a joule or more are necessary. Earlier experiments us- 
ing longer laser pulses also accelerated electrons via a 
self-modulation instability where the laser pulse length 
was much greater than the wavelength of a relativistic 
plasma wave, r > 27Tc/up [6HTQ]. However, these experi- 
ments were limited to operation at a low repetition rate 
(< 1 Hz) due to the high laser pulse energies. 

Using a low energy laser pulse to accelerate electrons 
requires a short underdense plasma, as laser depletion 
in generating the plasma wave limits the acceleration 
length. Conditions for electron trapping in the plasma 



wave are also restrictive due to the lower achievable inten- 
sity. One method for trapping electrons is to use a density 
downramp injection mechanism [11 , which was recently 
demonstrated experimentally using 10 TW lasers [12] [E]. 
In this scheme, the inhomogenous plasma leads to a time 
varying plasma wave phase velocity, which allows trap- 
ping once the electron velocity Ve exceeds the wave phase 
velocity Vph, Ve > Vph- 

In this Letter, we report on electron acceleration in 
an unexplored regime of plasma wakefield driven by few- 
millijoule femtosecond laser pulses (sub-TW) at high rep- 
etition rate (0.5 kHz). The high repetition rate enables 
better statistics and higher average particle flux, which 
were not accessible in previous experiments. Collimated 
electron beams are produced with quasi-monoenergetic 
spectra up to approximately 150 keV by acceleration in 
slow (non-relativistic) plasma waves on the density down- 
ramp of a 100 /im scale gas target. Because of the rel- 
atively high charge (~ 10 fC) and potentially short tem- 
poral duration, such electron sources have the potential 
to be used for ultrafast electron diffraction (UED) appli- 
cations. 

In conventional UED, electrons from femtosecond laser 
pulse induced photoemission are accelerated by an ex- 
ternal electric field. It is a challenging problem to con- 
trol the amount of charge and the temporal resolution of 
the electron bunch. Laser acceleration eliminates exter- 
nal acceleration instruments and stabilizes the electron 
bunch charge and duration. Recently, laser-accelerated 
electrons from solid target interactions have been demon- 
strated to successfully produce diffraction patterns [l^- 
Additional energy selection using bending magnets was 
crucial in Ref. [14 for obtaining mono-energetic electron 
pulses because the laser-accelerated electrons from solid 
targets have a very broad energy spectrum. In our exper- 
iments, the observed quasi-monoenergetic feature at 100 
keV can be exploited to allow a larger portion of electrons 
to be available for applications in electron diffraction and 
microscopy. The gas target also permits operation at 
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higher repetition rate with easier ahgnment. 

The experiments were performed using the laser 
system at the Center for Ultrafast Optical Science of the 
University of Michigan. This Ti: Sapphire based chirped 
pulse amplification (CPA) laser has a regenerative ampli- 
fier and a two-pass amplifier. It is capable of delivering 
pulses with energies up to 10 mJ and durations of 30 fs at 
a central wavelength of 800 nm. The output laser pulse 
is reflected from a deformable mirror (DM) and focused 
by either an //2 or //3 off- axis paraboloidal mirror to an 
optimized spot. Up to 8 mJ of pulse energy is available 
on target, which produces a peak intensity of 3 x 10^^ 
W/cm^ in a 2.5 /im (FWHM) focal spot. The laser focus 
is optimized by iteratively setting the DM so that the sig- 
nal of second-harmonic generation from a BBO crystal is 
maximized. 

To achieve high repetition rates with the desired small 
diameter gas profile, we used a free flowing capillary 
source. The gas target was produced by flowing pure 
argon or helium with a 5% impurity of nitrogen through 
the fused silica capillary tubing having an inner diameter 
of 100 jam. An approximately 10 cm length of this tubing 
was connected to a standard compressed gas system. A 
motorized XYZ stage was used to manipulate the capil- 
lary tubing to an accuracy of 2 /im. The gas flow experi- 
ences free expansion into vacuum and different densities 
were achieved by varying the backing pressure. A Mach- 
Zehnder interferometer configuration, using a beam con- 
taining 2% energy split from the main beam and prob- 
ing the gas flow transversely, was used to measure the 
plasma density profile. The interferogram data provide 
diagnostics on both the plasma spatial distribution and 
its temporal evolution, by varying the probe delay time. 
A 2D electron density map was reconstructed via Abel 
inversion of the phase-shift data as shown in Fig. [T] 



ing a maximum magnetic field strength of 0.02 T at the 
midplane was used for this purpose. A 0.5 mm aluminum 
slit was placed vertically before the magnets. The elec- 
tron energy distribution was also measured with a cus- 
tom built magnet spectrometer [T5| having an entrance 
aperture with a diameter of 3 mm and a solid angle ac- 
ceptance of 1 msr. The electron spectra were recorded 
on image plates (FUJI BAS-SR 2025, calibrated in ^). 
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FIG. 1. (a) Lineout of electron on- axis density, (b) Recon- 
structed electron density map. 



FIG. 2. (color online), (a)-(b): Images of electron beam pro- 
files measured on Lanex screen (He:N2 mixture gas. Exposure 
time 0.5 s) (a) a typical collimated electron beam. The white 
solid curves are vertical and horizontal lineouts. The dotted 
contours (half-maximum) show the maximum fluctuation over 
six data acquisition sequences; (b) filamented electron beam. 

(c) -(d): Lanex images showing the electrons deflected by the 
insertable magnets using He:N2 mixture (c) without magnets 

(d) with magnets, (e) Image plate raw data obtained using 
argon with the electron spectrometer accumulated over 2000 
shots. 



To measure the laser-accelerated electron beam profile, 
a Lanex scintillating screen was placed approximately 8 
cm from the nozzle along the laser propagation axis. The 
Lanex screen was covered by a 17 ±2 /im thick aluminum 
foil to block the laser light and cut off electrons with 
energies below ~ 50keV. Electron energies were mea- 
sured by horizontal magnetic deflection from their origi- 
nal propagation path and imaging on the Lanex screen. 
An insertable pair of disc magnets yoked together provid- 



Typical electron beam profiles are shown in Fig. [2|a) 
and (b). Well-collimated electron beams with a diver- 
gence angle (FWHM) as small as 2° were observed at 
lower densities [Fig. [2|a)]. As the backing pressure was 
increased and the nozzle position was lowered, the elec- 
tron beams broke up into several filaments [Fig. |2|b)]. 
When the laser was focused at a larger distance above 
the nozzle, the plasma length became greater due to gas 
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expansion and the peak plasma density decreased, which 
was confirmed by plasma interferometry. The overall ef- 
fect could be a longer interaction length at similar plasma 
densities leading to electron beam filamentation, as was 
also observed in laser wakefield experiments on 30 ter- 
awatt laser systems [17 . In both cases in these exper- 
iments, the features of the beam profiles remain stable 
when the system operated at 0.5 kHz. We also used an 
on-axis silicon PIN diode detector covered by aluminum 
foil to measure the shot-to-shot stability of the electron 
signal. The oscilloscope trace from the diode consistently 
showed less than 10% fluctuation. 

The raw images of the spectrum measurement are 
shown in Fig. |2|c)-(e) with calibration lines. The mea- 
sured spectrum exhibits a quasi-monoenergetic peak in 
the 100 keV energy range with width A£^fwhm ^ 20 
keV. In the experiments, data were obtained in "real 
time" for optimization of the beam parameters. For a 
fixed nozzle position, electrons were observed over a fi- 
nite range of backing pressures [Fig. [sj^a)] corresponding 
to plasma densities on the order of 10^^ cm~^, inferred 
from interferometric measurement (Fig.[T]). At lower den- 
sities, the plasma wave phase velocity is so high that the 
oscillating electrons are below the trapping threshold. At 
plasma densities where the plasma wake wavelength Xp is 
comparable to the laser pulse length L, large amplitude 
plasma waves are "resonantly" excited, enabling acceler- 
ation of fast electrons. 
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FIG. 3. (color online), (a) Measured electron beam charge 
and peak energy over a range of backing pressures, (b) Elec- 
tron charge measured by Lanex signal versus laser focus po- 
sition relative to the center of the capillary nozzle. Z < 
means laser is focused on the rear side of the nozzle. The 
blue line is the shape of a Gaussian density profile with a 
FWHM of 120 /xm. 



At higher densities, the laser pulse is susceptible to 
plasma defocusing or filamentation instability. The 
Lanex signal showed that electrons were preferably pro- 
duced when the laser was focused on the rear side of 
the nozzle [see Fig. [sj^b)]. This is likely to be related to 
the acceleration mechanism based on density downramp 
injection [12 , which will be discussed later. In the exper- 
iments, a gas mixture of helium and nitrogen was used 
under the consideration that electron trapping may be 
enhanced by additional ionization of the inner shell elec- 



trons near the peak laser intensity from a higher atomic 
number gas (e.g. nitrogen) [18 . However, the standard 
ionization induced trapping is not observed in our exper- 
iments, although the gas mixture yields a higher electron 
density compared to a pure helium gas with the same gas 
profile. 

To study the acceleration mechanisms, we performed 
numerical simulations using the 2D particle-in-cell (PIC) 
code OSIRIS [19]. The simulation ran in a stationary 
window of 814 jam x 38 /im, with a grid size of 22000x600 
cells. A Gaussian profile of neutral helium gas was used 
with the peak centered at 200 /im [see the lineout in Fig. 
|4ja)]. The peak density was O.OOSnc, where the plasma 
critical density Uc is 1.7 x 10^^ cm~^ for 800 nm light. 
The peak width (FWHM) was 120 /im as determined 
from the interferometric measurement (Fig.jl]). Electrons 
were produced using the ADK ionization model [20 , with 
4 particles per cell in each dimension (i.e. 16 fold). The 
laser parameters were chosen to match our experiment 
and consisted of a Gaussian spatial profile with a waist 
of wq = 1.8 /im, and a 5th order polynomial temporal 
profile similar to a Gaussian with a pulse duration of 
tp = 32 fs. The laser pulse leading edge was initialized at 
25 /im and focused at 200 /im. The simulation ran for 3.4 
ps. 

The short laser pulse generates large amplitude plasma 
waves as it propagates through the center of the gas from 
its ponderomotive force, but not to wavebreaking ampli- 
tude. Some time after the laser pulse leaves the plasma, 
wavebreaking of the plasma waves is observed and elec- 
trons are trapped and accelerated. The reason for this 
trapping is that the plasma waves are formed on the 
downramp of the gaussian profile, and hence their phase 
velocity Vph is not constant in time. In a ID inhomoge- 
neous plasma, the wavenumber /c of a plasma wave varies 
in time according to dkp/dt = —duOp/dx [21 . For an ap- 
propriately directed traveling wave on a density ramp, 
kp increases so that Vph reduces. When the wave phase 
velocity falls below the maximum electron oscillation ve- 
locity in the wake, the charge sheets will cross and trap- 
ping (wavebreaking) will commence [IJ^. Electron trap- 
ping in a gradual density inhomogeneity can occur several 
plasma periods behind the laser pulse [23 . This effect is 
observed in the simulations and illustrated in Fig. [4ja) 
and (b) , where the Fourier transform of the electron den- 
sity along the central axis is plotted against time, along 
with the number of trapped electrons. The peak elec- 
tron density Ue = O.Olnc in this simulation translates to 
a relativist ic plasma wavenumber kp = 0.1 /c^ for a lin- 
ear plasma oscillation, where kL = 27tc/Xl is the laser 
wavenumber in vacuum. 

In Fig. [4]^b), we plot the electron phase space show- 
ing the electrons injected in the accelerating phase of the 
wake. At 270 /im where trapping occurs, the phase ve- 
locity of the plasma wake in the simulation is estimated 
to be 0.3c. Using the cold, non-relativistic upper limit 
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for the wave breaking in the one-dimensional approxi- 
mation Emax = ^e^pVph/e [22 , the Calculated value 
is eEmax/'^e^LC ^ 0.02, which is slightly larger than 
the simulated value 0.015 [see Figgb)]. This is due to 
an equivalent thermal effect in the 2D simulations - the 
plasma oscillation has a spread in the fluid speed as seen 
in Fig|4]^b). Equivalent ID PIC simulations we ran con- 
firm that the trapping condition agrees with the analyt- 
ical expression very well. 

These trapped electrons are subsequently accelerated 
to a sub-relativistic forward momentum in the slow wave. 
At later times, a longitudinal electric field established by 
the space charge separation pulls electrons with lower 
energies back to the plasma but the energetic portion 
of the trapped electrons escapes as shown in Fig. [4]^c). 
To compare the electron spectra in the simulations with 
the experiment, the energy distribution of all electrons 
traveling in a forward cone angle of 50 mrad is plotted 
in the same graph with the measured spectrum in Fig. 
[5] The simulation result shows good agreement within 
the parameter uncertainty. The gas densities in the sim- 
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FIG. 5. (color online). Measured and simulated electron 
spectra. Black line with crosses: measured spectrum from 
Fig. |2je). Dashed blue line: energy distribution of all elec- 
trons in a forward cone angle of 50 mrad obtained from 2D 
PIC simulations. 



ulations were varied over a factor of 4 and ~ 100 keV 
electrons were observed over the whole density range, 
with slightly higher energies at lower densities, consis- 
tent with the experiments. The simulations suggest that 
the electron bunch has a few tens of fs duration close to 
the source. 

In conclusion, we have used a high repetition rate 10 
mJ short-pulse laser to demonstrate plasma wakefield ac- 
celeration of electrons. Reproducible collimated electron 
beams with a quasi-monoenergetic spectrum in excess of 
100 keV can be produced, for our parameters. Simu- 
lations show the electrons are trapped and accelerated 
to sub-relativistic energies in slow plasma waves. With 
the capability of operating the system at 500 Hz, fur- 
ther optimization and control of the electron beam pro- 
file, energy and stability can be realized. In addition to 
demonstrating the scalability of wakefield acceleration to 
lower energies, such a source may be useful for ultrafast 
electron diffraction applications. 
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FIG. 4. (color online), (a) The evolution of plasma wave 
number spectrum. Solid curve (red) is the total charge of 
the trapped electrons defined as those with pe/rrieC > 0.4. 
(b) Electron p-x phase-space plot at t = 1340 fs showing the 
onset of electron trapping. Dashed blue curve is a lineout of 
the Gaussian density profile. Solid red curve is the on axis 
longitudinal electric field, (c) Electron p-x phase-space plot 
at t = 3390 fs. Solid curve (red) is the averaged longitudinal 
electric field in the central axis region. Note the colormaps 
are in logarithmic scale. 
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